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ABSTRACT 

We analyze the dependence of the stellar disc flatness on the galaxy morphological 
type using 2D decomposition of galaxies from the reliable subsample of the Edge-on 
Galaxies in SDSS (EGIS) catalogue. Combining these data with the retrieved mod¬ 
els of the edge-on galaxies from the Two Micron All Sky Survey (2MASS) and the 
Spitzer Survey of Stellar Structure in Galaxies (S^G) catalogue, we make the following 
conclusions: 

(1) The disc relative thickness Zq/H in the near- and mid-infrared passbands cor¬ 
relates weakly with morphological type and does not correlate with the bulge-to-total 
luminosity ratio B/T in all studied bands. 

(2) Applying an ID photometric profile analysis overestimates the disc thickness 
in galaxies with large bulges making an illusion of the relationship between the disc 
flattening and the ratio B/T. 

(3) In our sample the early-type disc galaxies (SO/a) have both flat and “puffed” 
discs. The early spirals and intermediate-type galaxies have a large scatter of the disc 
flatness, which can be caused by the presence of a bar: barred galaxies have thicker 
stellar discs, on average. On the other hand, the late-type spirals are mostly thin 
galaxies, whereas irregular galaxies have puffed stellar discs. 
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1 INTRODUCTION 

Morphological and photometric characteristics of galaxies 
vary along the Hubble sequence (see the review by Roberts 
& Haynes 1994). Contribution of the bulge to the total lumi¬ 
nosity of a galaxy decreases, galaxies tend to be bluer, spiral 
pattern gets less tightly wound, the fraction of gas increases 
while going from early- to late-type galaxies. 

Guthrie (1992) analyzed the distribution of isophotal 
axial ratios in a large sample of UGC galaxies to study true 
axial ratios b/a and found the dependence of the true ax¬ 
ial ratio on the morphological type and that the late type 
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galaxies have smaller b/a. Karachentsev et al. (1997) mea¬ 
sured the apparent axial ratios 6/a of galaxies from the Flat 
Galaxies Catalog (Karachentsev et al. 1993) and found that 
the flat galaxies exhibit the dependence of the flattening on 
the morphological type: galaxies of the late types are flatter. 

Although the apparent axial ratio is an indicator of the 
disc thickness, these two characteristics are not the same, de 
Grijs (1998) performed the ID photometric profile analysis 
for 45 edge-on galaxies and derived global galaxy parameters 
in three photometric bands (optical and near-infrared). He 
found a correlation between the ratio of the radial to vertical 
scale parameter and galaxy type: galaxies become systemat¬ 
ically thinner when going from the type SO to Sc. According 
to de Grijs (1998), the ratio of the radial scale length h to 
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the vertical scale height zq in his sample of edge-on galaxies 
varies from 1.5-2 in the early-type spirals to 3-8 in the Sc- 
Sd galaxies. Note, however, that the scatter in zo/h is large 
for any certain morphological type. 

Kregel et al. (2002) re-analyzed the I band photom¬ 
etry of 34 edge-on spiral galaxies from the sample by de 
Grijs (1998). They applied a 2D algorithm to derive the 
scale parameters of the stellar disc and, thereby, the disc 
flattening. They also used the revised Hubble types of the 
galaxies from the HyperLeda database, which were signif¬ 
icantly shifted around as compared with those used by de 
Grijs (1998), sometimes by more than one subtype. How¬ 
ever, Kregel et al. (2002) did not use the revised data to 
discuss the dependence of the ratio Z(,/h on the morpholog¬ 
ical type found earlier by de Grijs (1998) on the base of the 
same sample. It was done by Hernandez & Cervantes-Sodi 
(2006), who analyzed the data by Kregel et al. (2002) and 
confirmed a diminishing trend of the ratio Zo/h while going 
from the early to late types: the late-type galaxies at the 
Sc end have thinner discs than the earlier types, on average. 
Gareful examination of fig. 7 in Hernandez & Cervantes-Sodi 
(2006) shows that the trend is supported by only one point 
(a galaxy of Sa type). Therefore, we can conclude that the 
reported correlation is rather poor. 

Ma et al. (1997, 1998, 1999) explored some statistical 
correlations for spiral galaxies of intermediate inclinations. 
They measured the disc thickness by the method proposed 
by Peng (1988). The method for determining the thickness 
of a spiral galaxy is based on the solution of Poisson’s equa¬ 
tion for a logarithmic disturbance of density. Ma et al. (1997, 
1998, 1999) showed that the ratio of the vertical to radial 
scales in the stellar disc tends to decrease along the Hub¬ 
ble sequence. It is important to note that all mentioned re¬ 
sults are based on indirect methods for determining the disc 
flatness, on a simple ID photometric profile analysis or on 
limited data sets. 

Bizyaev et al. (2014) (BKM14 hereafter) created the 
catalogue EGIS, a catalogue of edge-on galaxies in the Sloan 
Digital Sky Survey (SDSS). In their 3D photometric analysis 
of 5747 edge-on galaxies selected from SDSS Data Release 8 
they did not hnd significant correlation between the disc flat¬ 
ness and the morphological type. Moreover, a visual inspec¬ 
tion of galaxy images shows that early-type galaxies may 
have very flattened, as well as “puffed” discs (see Fig. 1). 
Bizyaev & Kajsin (2004) also noticed the existence of very 
thin stellar disks in galaxies with large bulges. 

To sort out this question, we analyze the dependence 
of the stellar disc flatness on the galaxy morphological type 
using robust 2D decomposition results obtained for a rep¬ 
resentative subsample of the EGIS catalogue. We also use 
the results of studying edge-on galaxies from the 2MASS 
archive (Mosenkov, Sotnikova & Reshetnikov 2010, MSRIO 
hereafter) and from the Spitzer Survey of Stellar Structure 
in Galaxies (S"^G, Sheth et al. 2010). 

This paper is organized as follows. In Section 2, we de¬ 
scribe our methodology of the 2D bulge/disc decomposition. 
We show that the approach used for decomposing edge-on 
galaxies brings up statistically reliable results. In Section 3, 
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Figure 1. Examples of galaxies from the EGIS catalogue (i 
band): IC 3608, Sb (Sab as given in the Catalogue of Detailed Vi¬ 
sual Morphological Classifications for 14,034 Galaxies in SDSS, 
Nair & Abraham 2010), z^jh = 0.10, bulge-to-total luminosity 
ratio BjT = 0.77 (top); and NGC 1032, SO/a (Sa in the Upp¬ 
sala General Catalogue of Galaxies, Nilson 1973), zojh = 0.30, 
B/T = 0.70 (bottom). The morphological types are taken from 
the HyperLeda database. The parameters are estimated via the 
2D bulge/disc decomposition using DEC A package (Mosenkov 
2014) with masking the central dust lanes (see Section 3.1 and 
also Fig. 4). 

we introduce the samples of galaxies which is considered in 
this work. In Section 4, we present the results of our anal¬ 
ysis, namely the correlations between the disc flatness and 
other quantities such as the morphological type, bulge-to- 
total luminosity ratio, galaxy colour, and galaxy luminosity. 
In Section 5, we discuss why the results of this article con¬ 
tradict the results of previous investigations. We also make 
some comments on the physical sense of the lack of the corre¬ 
lation between the disc flattening and the galaxy morpholog¬ 
ical type. In Section 6, we summarize our main conclusions. 
The cosmological framework adopted throughout this paper 
is Ho = 70 kms“^ Mpc“^, Q,m = 0.3, and Da = 0.7. 


2 THE DECOMPOSITION TECHNIQUE 

We first describe the main decomposition technique em¬ 
ployed in this study to analyze the structure of edge-on 
galaxies. There are three methods for decomposing such 
objects widely used in the structural analysis: the one¬ 
dimensional (ID, e.g. de Grijs 1998; Bizyaev & Mitronova 
2002 and BKM14), two-dimensional (2D, see Simard et al. 
2002; Peng et al. 2010; de Souza, Gadotti & dos Anjos 
2004), and three-dimensional method (3D, see, for exam¬ 
ple, Xilouris et al. 1999; Bianchi 2007; Baes et al. 2010 and 
BKM14). 

The ID methods applied to edge-on galaxies considers 
photometric profiles along the major and minor axes inde¬ 
pendently. In contrast to that, in 2D fitting the whole infor¬ 
mation from the 2D image is used simultaneously to build a 
robust model for each galactic component. There are several 
examples in the literature showing that the 2D method is 
much more reliable than the ID method (e.g. de Jong 1996) 
in retrieving the accurate structural parameters. In the 3D 
method a dust extinction model is usually adopted, which 
makes this method more precise but rather time consuming. 
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Since proper treatment of radiative transfer models takes 
significant computational resources (Bianchi 2007; Baes et 
al. 2010; De Looze et al. 2012), various simplifications can 
be applied (BKM14). In this section, we briefly describe the 
adopted 2D technique in the comparison with the ID and 
3D methods in Section 5.1. 


2.1 The 2D decomposition 

The 2D photometric decomposition of galactic images is re¬ 
alized in a number of packages such as GIM2D (Simard et 
al. 2002), GALFIT (Peng 2002; Peng et al. 2010), BUDDA (de 
Souza, Gadotti & dos Anjos 2004), and IMFIT (Erwin 2014). 
In this work we use three of them: BUDDA, IMFIT, and the 
most famous code GALFIT. Although the minimization algo¬ 
rithms are different in these packages, we assume they bring 
similar results of the decomposition since they have been 
tested and compared by different authors (e.g. Haussler et 
al. 2007, 2013; Kim et al. 2014; Busch et al. 2014). Also, 
each of these codes has its advantages which will be noticed 
further. 

To perform the bulge/disc decomposition, new Python 
wrapper DECA^ was made, as described in detail by 
Mosenkov (2014). The DECA software package was designed 
to perform photometric analysis of the structure of regu¬ 
lar galaxies, with the simple automatic decomposition onto 
bulge and disc using GALFIT as the galactic model genera¬ 
tor. It also takes advantages of several widely used packages 
such as IRAF, SExtractor and several Python libraries as, for 
example, PyFITS. DECA requires minimal human interven¬ 
tion as no any initial values of the bulge or disc parameters 
should be given. The algorithm is built in the way that it 
uses the input image of the studied galaxy (or, it may be a 
held with many different objects), the Point Spread Function 
(PSF) image, and the input text hie with the main informa¬ 
tion about the object (equatorial coordinates RA and DEC, 
redshift, and colour) and the frame calibration (photometric 
zero-point, pixel scale, gain, etc). 

A simple model of an edge-on disc galaxy represents a 
superposition of two photometric components: the bulge and 
the edge-on disc. The distribution of the surface brightness 
in the radial r and vertical « directions for the transpar¬ 
ent ‘exponential’ disc observed at the edge-on orientation is 
described by the following expression: 

I{r,z) = I{0,0)^Ki (^0 sech^(2/2o), (1) 

where I (0, 0) is the disc central intensity, h is the radial scale 
length, zq is the ‘isothermal’ scale height of the disc (Spitzer 
1942; van der Kruit & Searle 1981a,b, 1982a,b), and Ki is 
the modihed Bessel function of the hrst order. To hnd the 
edge-on disc central surface brightness, we should calculate 
j-edge-on _ l(r,0)dr (designated as /io,d if expressed in 
mag/arcsec^). 


^ http://lacerta.astro.spbu.ru/?q=node/96 
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The bulge profile is given by the Sersic law (Sersic 1968): 

/(r) = 7o,be-‘'"(’-/’'-")'^", (2) 

where Fe.b is the effective radius, or the half-light radius, /o.b 
is the central surface brightness (it can be reduced to the 
effective surface brightness /te.b expressed in mag/arcsec^), 
n is the Sersic index defining the shape of the profile, and 
the parameter depending on n ensures that re,b is the 
half-light radius. 

In order to estimate the bulge-to-total luminosity ratio 
B/T, we can find the total luminosity of the bulge and the 
disc via integrating the equations (1) and (2) over the whole 
galaxy image. 

DECA also allows us to take into account the stellar disc 
truncation when needed. The axis ratio of the bulge qh and 
its ellipticity parameter that controls the shape of the bulge 
(“boxy” or “discy”) can be set as free parameters. 

2.2 Robustness of the DECA decomposition 

As the main aim of this work is studying edge-on galaxies 
selected from SDSS, we tested DECA on the sample of ar¬ 
tificially created galaxy images that imitate SDSS images 
in the i band with the typical for the EGIS galaxies signal- 
to-noise ratio. We generated a galaxy image as the sum of 
photometric disc and bulge components by using expres¬ 
sions (1) and (2), respectively. We built a sample of 1000 
artificial images with the uniform distributions of the param¬ 
eters similar to those from Kregel et al. (2002). The bound¬ 
ary values for each parameter are taken from the physical 
background summarized in several works (e.g. Gadotti et 
al. 2009; MSRIO). The disc parameters cover the following 
ranges: 

16 /io.d 21 mag arcsec”^ , 

0.5 h lOkpc , 

0.1 zo SS 1.5 kpc. 

For bulge the ranges are: 

18 /Ue.b ^ 22 mag arcsec”^ , 

0.1 Te.b < 2 kpc , 

0.5 n ^ 5 , 

0.5 s; (?b sf 1.0. 

The image scale varied randomly between 0.4 and 2 
kpc/arcsec. The simulated images of galaxies were created 
using the Python wrapper SIGAL^ which operates the GAL¬ 
FIT code to produce model images. All created images were 
convolved with the Moffat PSF with FWHM = 1.2 arcsec. 
The pixel size is 0.396 arcsec/pix, the exposure time was set 

^ https://github.com/latrop/SIGAL 
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to 53.9 sec and the gain is 4e“/ADU. We added gaussian 
readout noise and poison noise to each pixel of the image. 
All galaxy discs were truncated up to the distance of 4/i. 

The decomposition of the artificial galaxy images was 
then performed with the DECA package. The robustness of 
the decomposition is shown in Fig. 2 for the parameters 
B/T and zo/h which are of great interest for us in this pa¬ 
per. We can see that for galaxies with a semi-major axis 
(being estimated as the truncation radius of the galaxy) 
larger than 20 arcsec {{AB/T)/{B/T)) = 0.01 ± 0.23 and 
{{Azo/h)/{zo/h)) = 0.017±0.06 while for galaxies with the 
lesser semi-major axis {{AB/T)/{B/T)) — 0.077±0.55 and 
{{Azo/h)/{zo/h)) = 0.05±0.18. Therefore we conclude that 
the minimum diameter of a galaxy with B/T > 0 to be 
well-decomposed is Dmin ~ 40 arcsec (hereafter we put the 
upper index to designate the passband). For galaxies with 
the diameter larger than Umin we may expect that the av¬ 
erage errors for the fitted parameters zo/h and B/T should 
be less than 10% and 25% of their true values respectively. 
Thus, the results of our test conhrm the robustness of fit¬ 
ting performing by DECA and its possibility for decomposing 
large sets of galaxies imaged by SDSS. 


3 THE SAMPLES AND THE 

DECOMPOSITION 

Our further analysis is based on the study of several samples 
of galaxies viewed almost perfectly edge-on (i ^ 86°) that 
were observed in different passbands. Although we focus our 
study on the sample of galaxies selected from SDSS, we also 
use the results of several investigations in near- and mid- 
infrared bands. 

The Sloan Digital Sky Survey (Aihara et al. 2011) is an 
ambitious project to investigate the Universe in five opti¬ 
cal passbands. Because the dust extinction along the major 
axis of the disc is substantially high in these bands, the 
vast majority of edge-on galaxy images demonstrate very 
well seen dust lanes. Hence, the proper model fitting of such 
galaxies should be carried out with taking into account the 
radiative transfer approach (e.g. Xilouris et al. 1999) which 
can be done for galaxies with sufficient spatial resolution 
(for nearby edge-on galaxies if we use SDSS) and takes sig- 
nihcant computational time. Nevertheless, as we will show 
later, some galaxies (of both early and late types) do not 
demonstrate sharp dust lanes in optics, and thus they can 
be decomposed using standard fitting packages described in 
Section 2. Galaxies with dust traces (having no sharp dust 
lanes) can be also decomposed using masking of these re¬ 
gions of high dust extinction. 

Since in the near-infrared (JHK) and mid-infrared 
Spitzer bands W1 (3.6 fim) and W2 (4.5 /rm) the flux is 
largely dominated by the stellar radiation rather than gas 
and dust emission (Whaley et al. 2009), and the dust extinc¬ 
tion here is much less than in optics, it is crucial to inves¬ 
tigate the structure of the stellar discs and bulges in these 
bands. Because of these aspects, the use of 2MASS (Skrut- 




Figure 2. Results of simulations as a function of the semi-major 
axis size of the artificially created edge-on galaxy images for the 
SDSS, i band, (a) The distribution of (B/T - B/T^t)/{B/T), 
where B/T is a known and B/Tf^^^ is a fitted bulge-to-total lu¬ 
minosity ratio, (b) The distribution of {zg/h — zq/ hf^^)/{ zq/K), 
where ZQ/h\s a, known and is a fitted disc flattening. The 

limiting radius is estimated to be about 20 arcsec (designated as 
the vertical red dashed line), or = 40 arcsec. The ordinates 

for each plot were averaged inside the bin of 5 arcsec (blue circles 
with the bars representing the standard deviation of the parame¬ 
ters inside each bin). The horizontal red dashed line corresponds 
to the zero value of the ordinate. 


skie et al. 2006) and S^G (Sheth et al. 2010) is suitable for 
our aims. 

In order to select edge-on galaxies from existing cata¬ 
logues and databases, we will use a selection criterion based 
on the galaxy apparent diameter at a given isophotal level. 
This allows us to select galaxies in the same manner. As 
was noticed in Section 2.2, the value of the limiting appar¬ 
ent galaxy diameter for each sample should be chosen so as 
to exclude galaxies with poor bulge/disc decomposition. In 
Section 2.2 we found that = 40 arcsec if we consider 

galaxy images in SDSS. Supposing that the main factor that 
affects the quality of bulge/disc decomposition is the image 
resolution and PSF (as reflected by the pixel size in arcsec). 
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we can estimate the minimum diameter of 2MASS and S'*G 
galaxies to be well decomposed by scaling our estimates per¬ 
formed for SDSS case. The pixel size in 2MASS is 2.5 times 
larger than in SDSS, hence the limiting apparent diameter 
should be ~ 100 arcsec. In the case of S^G where the 
pixel size is 0.75 pixel/arcsec, = 75 arcsec. 

In our study we consider only two main galactic compo¬ 
nents: the bulge and the disc, for simplicity. Possible influ¬ 
ence of the multi-component structure of edge-on galaxies on 
the main results of this work will be discussed in Section 5.3. 

The galaxy morphological type {Type) and the colour 
B — I for each galaxy were taken from the HyperLeda 
database (Paturel et al. 2003). The B and I bands are close 
to the g and i bands, respectively, therefore we opt to use 
the colour B — I for 2MASS and S'^G galaxies along with 
the colour g — i for SDSS galaxies. 

Each galaxy from the samples was visually inspected 
whether its shape resembles the shape of edge-on galaxies. 
Images in different passbands available in the NASA/IPAC 
Extragalactic Database (NED) were analyzed to check the 
presence of a dust lane along the edge-on disc. For all galax¬ 
ies of the samples we retrieved inclination angle i from the 
HyperLeda and found the mean value of all galaxies stud¬ 
ied (i) = 88.°0 ± 3.°5. In Appendix A we describe how we 
estimated the inclination angle in an independent way. 

Fig. 3 demonstrates the distribution of the total galaxy 
luminosities in our samples. We ensure that the galaxies 
used for the analysis are non-dwarf and luminous. Detailed 
description of the samples is given below. 

3.1 The SDSS sample 

In order to select comparatively large edge-on galaxies from 
SDSS, we used the EGIS catalogue which is described in de¬ 
tail by BKM14. This is the largest catalogue of verified edge- 
on galaxies at the moment which comprises 5747 genuine 
edge-on galaxies pre-selected automatically from the SDSS 
Data Release 8 and then visually inspected. This catalogue 
is complete for all galaxies with the major axis D* larger 
than 28 arcsec (see BKM14) and comprises disc galaxies of 
all morphological types (Sa-Sd-Irr), which is important for 
studying galaxies of different bulge contribution. 

In Section 2.2 we concluded that the minimum diameter 
of galaxies for keeping our decomposition robust is 40 arcsec. 
Therefore, from the EGIS catalogue we selected objects with 
D* > 40 arcsec and excluded galaxies with irregular struc¬ 
ture of discs and with various contaminants as bright over¬ 
lapping stars or galaxies. This resulted in the initial sample 
of 2027 objects. After a thorough visual inspection and pre¬ 
liminary decomposition of each galaxy image in the i band, 
we selected 145 galaxies and then estimated their inclina¬ 
tion angle using two methods described in Appendix A. We 
applied the hrst method to the galaxies without dust traces 
and the second method to the galaxies with visible dust man¬ 
ifestations. After that, we removed 4 galaxies with i < 86°. 
Therefore, the final sample includes 82 galaxies without any 
presence of dust lanes and 59 galaxies with slightly visible 
dust lanes. The mean inclination angle for the galaxies in 


our sample is (i) = 88.°2. An evidence that a galaxy has no 
strong dust attenuation came from the inspection of residual 
images (‘galaxy’—‘model’) and browsing of the vertical and 
radial photometric profiles. Since no dust traces were found 
in their profiles, we may hope that structural decomposition 
of these galaxies based on the model of transparent disc is 
more robust than for the rest of the sample. We will call this 
subsample of 82 ‘dust-free’ galaxies the “reference subsam¬ 
ple”. The selected sample consists of 27 early-type (T ^ 0) 
and 114 late-type galaxies (T > 0). The redshifts retrieved 
from the SDSS Data Release 9 prove that the galaxies of the 
sample are nearby {z < 0.05, (z) = 0.02). 

The 2D decomposition of the galaxy images from the 
sample of all 141 objects was performed using the DECA 
code. The i band images were downloaded from the SDSS 
Data Release 9 (Ahn et al. 2012). PSF (prefix “psField”) 
images were retrieved from the SDSS Data Archive Server 
(DAS). If a prominent dust lane was seen in the galaxy im¬ 
age (as for NGG 1032 in Fig. 4), it was masked during the 
decomposition. 

We checked our sample for completeness and found that 
the subsample of 98 galaxies with D* > 48 arcsec is com¬ 
plete while the whole sample is incomplete according to the 
E/Emax test (Thuan & Seitzer 1979). In Fig. 5, left panel, 
we can see that the distribution of the apparent axis ra¬ 
tio in the sample (gray color) and in the whole EGIS cat¬ 
alogue (red solid line) are very similar. We performed the 
Kolmogorov-Smirnov test to verify the null hypothesis that 
these two samples are drawn from the same population. We 
found that we cannot reject the null hypothesis at a 5% or 
lower a since the p-value is high and equals to 0.88. If the 
distributions of galaxy luminosity (Fig. 5, middle plot) and 
colour are considered, (Fig. 5, right plot) we can also see 
that these distributions look similar for both samples. How¬ 
ever, according to the Kolmogorov-Smirnov test there is no 
sufficient evidence to reject the null hypothesis (for the 
distribution: a = 0.105, p-value= 0.12; for the g — i distribu¬ 
tion: a = 0.09, p-value= 0.14). Since the apparent axis ratio, 
galaxy total luminosity and colour are the main photometri- 
cal characteristics of the observed objects, the distributions 
of these parameters for the sample of 141 galaxies and the 
EGIS catalogue are similar, we may conclude that our sam¬ 
ple well represents the whole EGIS catalogue. 

3.2 The 2MASS sample 

In order to create another sample of genuine edge-on galax¬ 
ies, we used a sample described in detail in MSRIO. The 
2MASS sample from MSRIO contains galaxies with a wide 
range of bulge contribution: from galaxies with prominent 
bulges to bulgeless galaxies. To select the galaxies, we 
used the 2MASS-selected Flat Galaxy Catalog (2MFGC, 
Mitronova et al. 2003) and the Revised Flat Galaxy Cat¬ 
alog (RFGC) by Karachentsev et al. (1999). Galaxies were 
selected from the 2MFGC according to the following cri¬ 
teria: the axial ratio sba > 0.2, the band Kron radius 

^ 30 arcsec and the concentration index (which is con¬ 
nected to the bulge-to-disc luminosity ratio) IC > 2.0. The 
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Figure 3. The distributions of the total luminosity for the galaxies from SDSS (i band), 2MASS (it's band) and S^G (3.6 /rm) samples. 



Model 
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Figure 4. The decomposition of the SO/a edge-on galaxy NGC 1032 using the DECA code (SDSS, i band). Left panel shows the 
galaxy image (top), the model (middle) and the residual image (‘galaxy’ —‘model’, bottom). Right panel represents the isophote map 
of the galaxy down to the 22.5 mag arcsec"^ surface brightness (top), the photometric cut along the major axis (bottom left) and the 
photometric cut along the minor axis (bottom right). The black dots refer to the galaxy, the red solid line represents the model, the blue 
dashed line represents the bulge, and the disc is shown by the green dots. The dust lane was masked during the decomposition. Main 
derived parameters of the decomposition are B/T = 0.70, zojh = 0.30. 


interacting, peculiar and non-edge-on galaxies identified vi¬ 
sually and according to the information found in HyperLeda, 
were removed from our consideration. The final sample con¬ 
sists of 175 galaxies in the Ks band, with images taken from 
the 2MASS archive. All the galaxies were decomposed into 
a bulge and a disc using the BUDDA code (de Souza, Gadotti 
& dos Anjos 2004) (see several examples of decomposition 
in MSRIO). 

We found that about one-third out of the sample of 175 
galaxies have inclinations less than 86° as estimated using 
the first method described in Appendix A. We removed these 
galaxies, as well as the galaxies with the diameter smaller 
than = 100 arcsec. The resulting sample comprises 66 
large genuine edge-on galaxies. 


According to the H/Hmax test our final sample is in¬ 
complete and thus may suffer selection effects. In the left 
panel in Fig. 6 we can see that the distribution of the fi¬ 
nal sample (gray color) and the distribution of the whole 
2MASS sample over the apparent axis ratio differ since we 
excluded non edge-on galaxies from subsequent analysis - 
these galaxies form a peak at b/a = 0.32 (the red line dis¬ 
tribution). The other distributions (middle and right pan¬ 
els) do not show significant differences between the initial 
sample of 175 galaxies and the final one. The selected sam¬ 
ple comprises 10 early-type and 56 late-type galaxies. All 
the galaxies are nearby, with 2 : < 0.035 (average redshift 
( 2 ) = 0.007); the mean inclination angle is {i) = 88.°0. 
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Figure 5. The distributions of the apparent axis ratio, total luminosity and colour for the selected SDSS sample (gray) and the EGIS 
catalogue (red solid line). 




Figure 6. The distributions of the apparent axis ratio, total luminosity and colour for the selected 2MASS sample (gray) and the sample 
from MSRIO (red solid line). 


3.3 The S^G sample 

The catalogue consists of 2352 nearby (with distance 

less than 40 Mpc) galaxies at 3.6 and 4.5 ^m, and it is 
a valuable source of edge-on galaxies. Using the Gal- 
ht Models Home Page®, we retrieved the best decomposi¬ 
tion models with the edge-on disc included in 362 galaxies 
(at 3.6/im). However, the visual inspection of the Spitzer, 
DSS and SDSS images of these galaxies showed that some 
galaxies are not quite edge-on since they have round ellip¬ 
tical isophotes, or visible spiral arms, or rings), therefore 
we did not classify these galaxies as edge-on (for example, 
NGG 5984 and NGC 678). Otherwise, we would risk to re¬ 
trieve false thicker stellar discs, which can signihcantly mis¬ 
lead us in studying the dependence between the disc flatten¬ 
ing and the galaxy morphological type. 

From the pre-selected sample of edge-on galaxies in 

® The catalogue of galactic structural parameters can be found at 
http:// WWW. oulu.fi/astronomy/S4 G-PIPELINE4 /s4 g-p4-table8. dat 


the S"^G model catalogue we selected galaxies consisted ei¬ 
ther of the disc alone (late types), or of the disc and the 
bulge (early and intermediate types'^). Galaxies with de¬ 
tected double-exponential discs are not considered in this 
work. Then, we estimated inclinations of the selected galax¬ 
ies using the hrst method from Appendix A and removed ob¬ 
jects with i < 86°. Also, we removed non-bulgeless galaxies 
with smaller than minimum diameter at the outer isophote 
of W1 = 25.5 magarcsec”® where was 

estimated to be 75 arcsec. 

The total number of edge-on galaxies of all types se¬ 
lected from the S^G model catalogue is 99, with the mean 
inclination angle (i) = 87.°9. The sample mostly consists of 
the late-type galaxies (8 early-type and 91 late-type galax¬ 
ies). We did not reject galaxies of the latest morphological 
types (23 irregular galaxies, with poorly defined structure) 
though we assume that their edge-on orientation determina- 


^ 31 galaxies were fitted with the ‘bulge-|-disc’ model 
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Figure 8. Examples of early-type galaxies from the S^G database 
(W1 band): NGC 1596 (SO, z^/h = 0.34, B/T = 0.64, top) and 
NGC 7814 (Sab, z^/h = 0.15, B/T = 0.90, bottom). The disc 
thicknesses and bulge-to-total luminosity ratios are taken from 
S^G Galfit Models Home Page. The red contours refer to the 
outer isophote of Wl=25.5 magarcsec“^. 

tion is rather precarious. We make some comments on this 
issue in the further sections. 

In Fig. 7 we show the distributions for the selected sam¬ 
ple (gray color) and the whole sample of 362 galaxies decom¬ 
posed with edge-on discs in the catalogue (the red line 
distribution). It is well seen that all present distributions are 
very similar. However, the distribution of shows that 
our S^G sample comprises larger number of low-luminosity 
galaxies and less number of high luminosity galaxies in the 
comparison with the whole sample of edge-on galaxies in the 
S'^G catalogue. 

Although our S'*G sample is significantly incomplete ac¬ 
cording to the H/Vlnax test and different morphological types 
are not statistically representative in the sample, this is the 
only sample of genuine edge-on galaxies observed in the mid- 
infrared bands (see examples of the edge-on galaxies from 
the created S^G sample in Fig. 8), thus we consider this 
sample carefully in the current work. 


4 RESULTS 

4.1 The disc thickness versus the morphological 
type 

First, let us consider the dependence of the disc relative 
thickness zo/h on the galaxy type. In hg. 14 from BKM14 
one can see the dependence between the stellar disc thickness 
and the type derived from the ID analysis. This dependence 
is similar to that found in de Grijs (1998) also based on the 
analysis of ID photometric prohles. At the same time, the 
3D analysis for the same sample (see fig. 15 in BKM14) does 
not allow us to conclude that this dependence truly exists. It 


only shows that the late-type edge-on galaxies (determined 
as galaxies without bulge) are thinner than others, on aver¬ 
age. Thus, we can conclude that a more advanced 3D method 
of decomposition erases the thickness-type dependence. 

It is worth noticing that morphological types in BKM14 
were estimated automatically for each galaxy using main 
photometrical parameters from SDSS. We compared these 
types with those extracted from the HyperLeda database 
and found the mean scatter between these two systems of 
morphological classification being 2 T-units. 

Fig. 9 shows the distribution of the galaxies from our 
SDSS, 2MASS and S'^G samples in the same zo/h-Type 
space. The SDSS sample is shown in Fig. 9, the top panel. 
We also marked here the galaxies with dust traces (open 
circles) and with visible X-shaped structure (red crossed 
symbols) which are candidates to hosting bars as discussed 
by Bureau et al. (2006). The whole sample demonstrates a 
trend (the Pearson’s correlation coefficient p = —0.56, the 
two-sided p-value p = 1.44e-12). The difference in the thick¬ 
ness between the thickest and the thinnest discs is of the 
factor of 1.5. 

For the 2MASS sample of large galaxies we do not see 
strong variations of the relative thickness with the morpho¬ 
logical type (p = —0.41, p = 0.001; see Fig. 9, middle plot). 
Since the S'^G sample consists mainly of late-type edge-on 
galaxies (many of them occur to be irregular), we have even 
an upward trend for the galaxies with T > 6 (blue dashed re¬ 
gression line with p = 0.48 and p = 0.016 on Fig. 9, bottom 
panel). For the galaxies with T ^ 6, the correlation between 
zo/h and Type is not seen (p = —0.15, p = 0.443). Notwith¬ 
standing small number of early-type galaxies and excluding 
irregular galaxies with T ^ 8, the overall picture does not 
reveal any correlation between the morphological type and 
the stellar disc flattening. 

The poor dependence between the disc flattening and 
the galactic morphological type is confirmed indirectly by 
Chudakova & Silchenko (2014). The latter paper presents a 
new photometric method for deriving the relative thickness 
of the stellar disc from 2D surface brightness photometry. 
The method was applied to a sample of 45 early-type (SO- 
Sb) galaxies. They found that the discs of lenticular and 
early-type spiral galaxies have similar thickness, on average. 

4.2 The disc thickness versus the bulge-to-total 
luminosity ratio 

Determination of the morphological type for edge-on galax¬ 
ies is difficult, van den Bergh (1998) in his review on the 
galactic morphology and classification stressed that edge-on 
galaxies are often difficult or impossible to classify. Naim et 
al. (1995) concluded that the agreement between indepen¬ 
dent observers to do morphological classihcation is good, on 
average, but the scatter is large (about 1.8 subtypes). How¬ 
ever, in the case of edge-on galaxies the agreement between 
independent classifiers is worse (see fig. 6, right panel in that 
work). In Section 5.2 we will demonstrate how the inclina¬ 
tion of a galaxy can affect its morphological classihcation. 

For edge-on galaxies it is more reliable to use the bulge- 
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Figure 7. The distributions of the apparent axis ratio, total luminosity and colour for the selected S^G sample (gray) and the whole 
sample of 362 edge-on galaxies in the S^G catalogue (red solid line). 


to-total luminosity ratio B/T as the indicator of the galaxy 
morphological type because this value can be estimated di¬ 
rectly from the photometric analysis and modeling. 

In Fig. 10 we plotted the dependence of the disc rel¬ 
ative thickness zo/h on the bulge-to-total luminosity ratio 
B/T. One can see that for galaxies from SDSS (the top 
panel) with 0.2 < B/T < 0.6 there is a large scatter of 
zo/h values around zo/h ~ 0.3. The galaxies with visible 
X-shaped structnres (red crosses) lie within the same lim¬ 
its 0.2 < B/T < 0.6 with the relative thickness zo/h ~ 
0.30 ± 0.05 in comparison with the average flatness for the 
whole sample (zo/h) = 0.25 ± 0.07. Thus, we may conclude 
that the galaxies with X-shaped structures have thicker stel¬ 
lar discs, on average. 

We verified whether our fitting procedure could af¬ 
fect the possibly existed correlation between zo/h and B/T 
taking into account the fact that these retrieved parame¬ 
ters may have 10% and 25% uncertainties respectively (see 
Sect. 2.2). We created a simulated sample with the same dis¬ 
tribution of zo/h as for the SDSS sample, where there exists 
correlation between the disc flatness and the bulge-to-total 
luminosity ratio with the regression coefficient p = 0.55. 
We then applied the biasing to each point in the B/T-zo/h 
space using the normal distributions with the standard de¬ 
viations 0.25 • B/T and 0.1 • zo/h for the values of B/T and 
Zo/h respectively. We repeated this biasing 100000 times and 
found that the probability to receive the final sample with 
the uncorrelated (p < 0.15) parameters zo/h and B/T is 
5 • 10“®. This test proves that the lack of the correlation 
between the disc flatness and the bulge-to-total luminosity 
ratio can not be caused by the fitting biases. 

We also marked by squares early-type galaxies with 
Type ^ 0 (Fig. 10, top plot). All these galaxies show a wide 
range of the ratio B/T from 0 to 0.8. They also demonstrate 
a significant scatter of zo/h whereas the ratio z^/h does not 
depend on B/T for these objects in accordance with the con¬ 
clusion by Chudakova & Silchenko (2014). The dependence 
in Fig. 9 (top plot) is mainly created by these galaxies. The 
classification of these galaxies as SO is expected to be biased 


because we can not conclude the absence of a spiral struc¬ 
ture in them. The number of galaxies classified as early-type 
galaxies in 2MASS and S^G samples, in turn, is small. As a 
result, the dependence of zo/h on Type for these samples is 
weak (Fig. 9, middle and bottom plots). 

For both 2MASS and samples we did not find any 
correlation between the flattening zq/H and the ratio B/T 
as well. From the shown distribution we can conclude again 
that the early-type galaxies may have stellar disc with dif¬ 
ferent flatness, including very thin discs. The correlation be¬ 
tween the relative thickness and the bulge contribution to 
the total luminosity is not confirmed for our samples. 


4.3 The disc thickness versus the galaxy colour 

One of the prominent correlations along the Hubble se¬ 
quence is the morphological dependence on the integrated 
colour (see the review by Roberts & Haynes 1994). The 
early-type galaxies are mostly red whereas the late-type ones 
are mostly blue. Thus, the colour, in addition to the ratio 
B/T, is another indicator of the morphological type. The 
dispersion of the relation between the colour and the mor¬ 
phological type is large: there are some galaxies of Sc type 
which are red, and vice versa, there may be blue SO/a galax¬ 
ies. As an example, among the sample of edge-on galaxies 
we found FGC 1597 (Scd) to be red and NGG 3600 (Sa) to 
be blue. 

In fig. 7 from BKM14 there is no significant trend in the 
disc thickness with the overall galaxy colour g — r. In Fig. 11 
we plot the z^/h versus the g — i for the SDSS sample (left 
plot). The colours were corrected for the Milky Way redden¬ 
ing using the extinction maps by Schlegel et al. (1998), but 
they were not corrected for the internal extinction in the 
galaxies. 

Fig. 11 (left plot) shows that there are two perpen¬ 
dicularly oriented clouds of galaxies. The bluer galaxies 
((/ — i < 1.0) have stellar discs thinner {zo/h ~ 0.21 ± 0.03) 
than the red galaxies, which additionally reveal a wide scat- 
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Figure 9. The distribution of the disc flatness as a function of 
the morphological type for the SDSS (top panel, i band), 2MASS 
(middle panel, it’s band) and S‘*G (bottom panel, 3.6 ^m) sample. 
The open circles in the top panel refer to the subsample of galaxies 
with slightly visible dust lanes (59 galaxies). The filled circles 
refer to the reference subsample of galaxies without dust lanes 
(82 galaxies). The red crosses mark the galaxies with sharp X- 
pattern structures in the central region. The blue lines represent 
the regression lines for each sample (see text). 


ter of the thickness {zo/h « 0.28±0.06). It is well seen that 
there are red galaxies with both flat and thickened discs. 

The latter can be interpreted in two ways. First, red 
galaxies (representing mostly early and intermediate types) 
can truly have thick or thin discs (as we could see in Fig. 1). 
Second, the dust lane (and the internal extinction in general) 
can cause the reddening. We tried to properly mask the dust 
lane while preparing the image for the decomposition. But 
remaining traces of the dust lane can alter the profile in the 
2 -direction making it shallower. It results in a larger value 
of the 2 o. In the latter case the results of the decomposi¬ 
tion should be biased, and the relative thickness zo/h should 
be slightly overestimated. Almost all galaxies with the dust 


Figure 10. The distribution of the stellar disc flatness as a func¬ 
tion of the bulge-to-total luminosity ratio for the SDSS {i band), 
2MASS {Kb band) and S^G (3.6 pm) samples. Symbols as in 
Fig. 9, and the green squares represent early-type galaxies with 
Type ^ 0. 

traces fall into the red cloud in the plot, while, as we can 
see, the dust-free galaxies fall into both clouds (Fig. 11, hlled 
circles). Note that we can not exclude biased decomposition 
in the i band even in the case of dust-free galaxies since we 
do not see the division into two clouds when using SDSS 
sample for the ratio B/T and Type. The galaxies with dust 
lanes sample the whole range of value B/T and Type. It 
indicates that the presence of the dust lanes marginally cor¬ 
relates with the B/T and Type (see, for example, Draine et 
al. 2007). 

For both 2MASS and S"^G samples there is no correla¬ 
tion between the disc thickness and the overall galaxy colour 
B — I (middle and right plots). We can see only a slightly 
larger scatter for the ratio zo/h in the redder side of the plot 
for the 2MASS sample. 

Thus, the difference between the distributions of the 
parameters estimated from the optical SDSS and IR sam- 
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pies on the zq / h-colour diagram may be interpreted by the 
presence of the dust component that affects the model de¬ 
composition. We suppose that even our ‘dust-free’ galaxies 
are not absolutely transparent in the i band, which biases 
the zo/h estimation. However, this issue should be consid¬ 
ered more carefully. 

Another possible factor that biases interpretations of 
the plots is that the S^G sample is deficient in early-type 
red galaxies and comprises a large fraction of blue galaxies. 
This fact can be one of the reasons why we do not see any 
subdivision of galaxies in the zo/h-B — I space for the S'’G 
sample. 

4.4 The disc thickness versus galactic luminosity 

We also analyze the dependence of the relative thickness 
zo/h on the total luminosity of galaxies and find this cor¬ 
relation is weak if at all existing (see Fig. 12). Our SDSS 
sample demonstrates a large scatter for the ratio zo/h in 
the bright side of the plot (Fig. 12, left panel), which is 
similar to the presence of two perpendicular clouds in the 
Zo/h — g — i diagram (see previous section). This scatter is 
mainly due to the presence of early-type galaxies {Type 0, 
squares in the plot), which show a large scatter of the ratio 

Zo/h. 

It is well established that the dust content correlates 
well with the galaxy mass (e.g. Masters et al. 2010; Skibba 
et al. 2011), and we see that the bright region is populated 
mainly by the galaxies with dust traces. 

Only one sample (S^G) shows a hint of the thickness - 
luminosity relation: more luminous galaxies are thinner, on 
average, de Grijs (1998) also did not find any evidence for a 
dependence between the relative thickness and the absolute 
magnitude (see fig. 7b in de Grijs 1998). He referred to the 
theoretical prediction that the ratio h/zo should decrease 
rapidly from faint galaxies to a constant level for normal and 
bright galaxies (Bottema 1993) and claimed that it could not 
be confirmed observationally in his sample. Our samples also 
do not confirm this trend. Moreover, S^G sample shows the 
inverse trend probably because of the large number of late- 
type (irregular) galaxies. However this trend is consistent 
with the issue of the thickening of disc towards low mass 
reported in recent years in the literature (Sanchez-Janssen 
et al. 2010; Roychowdhury et al. 2013). 

5 DISCUSSION 

5.1 The difference between the ID and 2D 
decomposition results 

The dependence of the stellar disc relative thickness on the 
Type reveals itself clearly only in the ID modeling (de Grijs 
1998; BKM14). This statement is also applicable to the de¬ 
pendence between the disc flatness zo/h and the bulge-to- 
total luminosity ratio B/T. Fig. 13 demonstrates this depen¬ 
dence for the SDSS sample in the case of ID analysis. The 
correlation between the disc thickness and B/T is good (the 
Pearson’s correlation coefficient p = 0.75) for the subsample 



Figure 13. The distribution of the disc flatness as a function 
of the bulge-to-total luminosity ratio for the whole EGIS sample 
(gray dots) and for the subsample of 141 galaxies (black circles) 
estimated using the ID decomposition in the i band. The dashed 
regression line is referred to the whole sample; the solid regression 
line is referred to the subsample. 


of 141 galaxies (black circles), and it is fairly good for the 
whole EGIS catalogue (gray dots, p = 0.53). However, the 
correlation does not appear at all if the 2D analysis results 
are used (see Fig. 10, top panel). 

The correlation between the zo/h and B/T parameters 
for the ID analysis of the SDSS subsample galaxies can be 
explained by biased estimation of the disc structural param¬ 
eters for the galaxies with significant bulge contribution. 
This is rather expected because in the early-type galaxies 
the extension of the bulge structure is much larger than 
that in late-type galaxies. The disc scale height zo should 
be overestimated while analyzing the vertical photometric 
disc profiles contaminated by the bulge light. At the same 
time, the radial scale length should be underestimated. As 
a result, the ratio zo/h will be overestimated in the early- 
type galaxies. The biased ID decomposition leads also to the 
underestimation of the ratio B/T in comparison with 2D de¬ 
composition. In Fig. 14 we compare the results for the zo/h 
and B/T approximation between the ID and 2D decomposi¬ 
tion techniques applied to our SDSS sample. The ratio zo/h 
is overestimated for the ID decomposition in comparison 
with the 2D decomposition (left panel). 

The biased decomposition in the ID analysis reveals it¬ 
self in the biased distribution of the ratio B/T (Fig. 13 and 
Fig. 14, right panel) in the comparison with the distribu¬ 
tion over B/T in the 2D analysis {B/T < 0.6 for the ID 
decomposition, and B/T spreads up to 0.8 for the 2D de¬ 
composition for the same galaxies). In the ID analysis the 
contribution of the bulge seems to be underestimated and 
the rest of the true bulge light (non-retrieved from the pro¬ 
file) contributes to the disc thickness. 

In Fig. 15 the ratio of the stellar disc thickness in the 
ID and 2D analysis versus the bulge-to-total luminosity ra¬ 
tio (from the 2D analysis) is shown. It is well seen that with 
the increasing of the B/T the relative thickness is growing 
with respect to the 2D approach (both for the whole sam- 
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Figure 11. The distribution of the disc flattening as a function of the colour g — i for the SDSS {i band), 2MASS {Ks band) and S'^G 
(3.6 fj.m) sample. The red and blue ovals outline the two possible families of galaxies selected from their loci on the g — i — zo/h plane. 
The symbol designation is kept the same as in Fig. 10. 





Figure 12. The dependence of the disc relative thickness on the 
(3.6 fim) samples. The designation is kept the same as in Fig. 10. 

pie and the reference subsample of 46 dust-free galaxies). 
In other words, the correlation in Fig. 13 may be a conse¬ 
quence of the biased model fitting in the ID analysis, which 
dramatically affects the decomposition results for galaxies 
with significant non-disc components. A simple exponential 
approximation of data is shown in Fig. 15 by dashed line. We 
can see that for galaxies with B/T > 0.4 the ID disc flat¬ 
tening will overestimate more than 10% the flattening found 
with the the 2D method and more than twice for galaxies 
with B/T ~ 0.8. 

We summarize that the ID and 2D methods produce 
different results as we can see by comparing Fig. 13 and 
Fig. 10 (top panel). We assume that this difference is caused 
by biased estimation of the bulge contribution to the to¬ 
tal luminosity in the ID approach. With the ID analysis, 
the bulge was not fitted at all and only its fraction was es¬ 
timated which may bring wrong decomposition results for 


total luminosity for the SDSS {i band), 2MASS (Ks band) and S'^G 


bulge-dominated galaxies. The 2D approach is a more ro¬ 
bust method, which should report reliable results. The 3D 
analysis performed in BKM14 seems to be rather ambigu¬ 
ous, and it should be properly improved. Therefore, we do 
not discuss the results of this method here. 

5.2 Morphological type estimation 

The morphological classification of edge-on galaxies is a dif¬ 
ficult task because of several aspects. First, the most essen¬ 
tial criterion of the Hubble type, the tightness of the spiral 
pattern, cannot be applied to galaxies viewed edge-on. Sec¬ 
ond, dust attenuation dims the bulge in optics and NIR and 
hence, we are not able to estimate by eye the contribution 
of the bulge to the total luminosity in the galaxy reliably. 
Third, some structural details in a galaxy can be identified 
only via the inspection of its non-edge-on surface brightness 
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Figure 14. The comparison of the results for the zg/h and B/T approximation between the ID and 2D decomposition techniques. The 
open circles refer to the subsample of galaxies with slightly visible dust lanes; the filled circles refer to the reference subsample of galaxies 
without dust lanes. 




Figure 15. The dependence of the ratio of the stellar thicknesses 
in the ID and 2D modeling versus the bulge-to-total luminosity 
ratio from the 2D analysis for our SDSS sample, i band. The open 
circles refer to the subsample of galaxies with slightly visible dust 
lanes; the filled circles refer to the reference subsample of galax¬ 
ies without dust lanes. The dashed line is a simple exponential 
approximation to the data. 

maps. For instance, steps (“shoulders”) at the photometric 
profile along the major axis of an edge-on galaxy can be an 
evidence of a bar or a ring, but they can also reveal the 
deviation of the disc surface brightness profile from the sin¬ 
gle exponent. Wrong morphological classification has been 
found for many galaxies when detailed bulge/disc decompo¬ 
sition was performed (see, for example, de Souza, Gadotti 
& dos Anjos 2004). 

In order to demonstrate how the inclination of a galaxy 
can affect its morphological classification, we show the de¬ 
pendence between the B/T and Type for two subsamples 
from the S'^G catalogue (Fig. 16). The first group contains 
galaxies with i < 30° (black circles). The second one com¬ 
prises galaxies with i = 90° (open circles). Here, the inclina¬ 


Figure 16. The distribution of the bulge-to-total luminosity ratio 
as a function of the morphological type for galaxies with inclina¬ 
tion angles i = 90° (open circles) and i < 30° (black circles), in 
the S^G sample (3.6 pm). 

tion angle i is taken from the HyperLeda database. Although 
some of these galaxies are not viewed perfectly edge-on, they 
are certainly highly inclined. In Fig. 16 one can see that 
the open circles have larger scatter, with many outliers in 
the right upper corner. This indicates that the morpholog¬ 
ical classification of edge-on galaxies is rather ambiguous. 
Observable strong dust attenuation can dim the bulge and 
hence bias the morphology classification for such objects to¬ 
wards the later types. 

5.3 The model dependence of the obtained results 

We should now make some comments on decomposing mul¬ 
ticomponent galaxies. It is obviously difficult to unam¬ 
biguously identify the presence of a bar in nearly edge- 
on systems. However, about half of disc galaxies demon- 


(c) 2015 RAS, MNRAS 000, 1-16 










14 A.V. Mosenkov, N.Ya. Sotnikova, V.P. Reshetnikov, D.V. Bizyaev, and S.J. Kautsch 


strates boxy/peanut-shaped structures in the central region 
(Liitticke et al. (2000); Laurikainen et al. (2011)), which is 
tagged by some authors as projection of the vertical struc¬ 
ture of bars (Bureau & Freeman 1999; Chung & Bureau 
2004). 

It is apparent that for such galaxies the simple 
‘bulge-|-disc’ model will not provide appropriate fitting. In 
this case a ‘bulge-l-disc-l-bar’ modeling should be used in¬ 
stead. The bulge-to-total luminosity ratio and the relative 
disc thickness are expected to decrease in the case of the bar 
contribution. More detailed images of nearby galaxies show 
some other hidden features. For instance, the photometric 
decomposition of the Sombrero galaxy NGC 4594 revealed 
a disc, a bulge, an outer spheroid, a stellar ring, and a stel¬ 
lar, inner ring or disc (Gadotti & Sanchez-Janssen 2012). 
However, in the case of such multi-component structure the 
fitted model becomes ambiguous and suffers from the local 
minimum degeneracy. 

Also, some authors tried to find the parameters of thin 
and thick discs using a two-component model of the stel¬ 
lar disc (see, e.g. Yoachim & Dalcanton 2006; Gomeron et 
al. 2011). In this work we only studied one-component disc 
decompositions (mainly because of the shallow photometric 
depth of the surveys we work with and a lack of sufficient 
spatial resolution) following the works by de Grijs (1998) 
and Kregel et al. (2002). Studying larger samples in different 
photometrical bands, we did not find a correlation between 
the disc flattening and the morphological type, as well as 
between the disc flattening and the bulge contribution to 
the total galaxy luminosity. 

Our next step in studying edge-on galaxies is building 
the best models of nearby, well-resolved edge-on galaxies us¬ 
ing their multi-band imaging, with proper dust-component 
treatment (see e.g. Baes et al. 2011). 


6 CONCLUSIONS 

The results of our analysis of edge-on galaxies can be sum¬ 
marized as follows: 

(i) In the course of our analysis we obtained a counter¬ 
intuitive result, — the disc relative thickness does not de¬ 
pend on the bulge-to-total luminosity ratio for all the sam¬ 
ples studied. However, it may depend on the total luminosity 
in the near- and mid-infrared bands, in which the old stellar 
population contributes mainly in the disc luminosity, and 
reflects the dependence on the total mass (Sanchez-Janssen 
et al. 2010; Roychowdhury et al. 2013). The existence of the 
correlation between the disc flatness and the morphological 
type is only confirmed for the SDSS sample and late-type 
galaxies from the S‘*G sample (anti-correlation). 

These findings may have far-reaching consequences. Thus, 
since the relative thickness of the disc depends on the con¬ 
tribution of the spherical component (‘bulge-fdark halo’) 
{h/zo oc Mdark/Mtotai; Zasov, Makurov & Mikhailova 1991; 
Zasov et al. 2002, see also the discussion in MSRIO and 
Mosenkov, Sotnikova & Reshetnikov 2014), this means that 
the relative contribution of the spherical component does 


not vary systematically along the morphological Hubble se¬ 
quence of bright spiral galaxies. In turn, since bulges con¬ 
tribute more to the mass of the spherical component in the 
early-type disc galaxies, the dark halos themselves should re¬ 
veal higher mass fraction in the case of the late-type (mostly 
bulgeless) galaxies. 

At the same time, a compact but massive bulge stops the 
disc heating in the vertical direction due to the bending 
instability (Sotnikova & Rodionov 2006) and leaves the disc 
rather thin. In galaxies with large B/T ratio we may expect 
observing stellar discs of either thickness. In Fig. 10 we see 
a large scatter of the ratio zo/h for the galaxies with B/T > 
0.4. 

(ii) We conclude that the correlation between the mor¬ 
phological type and the disc flatness reported earlier by some 
authors (e.g. de Grijs 1998) was due to biased estimates of 
the disc flatness during the ID analysis. Reanalyzing the 
same data by the 2D method Kregel et al. (2002) found the 
disc scale height to be overestimated for almost all galaxies 
studied with the ID decomposition (see their fig. 2, bottom 
panel). At the same time, the disc scale length was underes¬ 
timated (the same fig. 2, top panel). As a result the relative 
thickness increases in the the ID analysis. We obtained the 
same result and extended it: the ID decomposition overes¬ 
timates the disc thickness mainly in the galaxies with large 
bulges (more than 10% for galaxies with B/T > 0.4), i.e. 
in all early-type disc galaxies. It is important to note that 
Kregel et al. (2002) did not present the correlation between 
the disc thickness and the galaxy type for the updated anal¬ 
ysis. It was done by Hernandez & Cervantes-Sodi (2006) 
and is shown in their fig. 7. A glance at the plot suggests 
that the reported trend is poor and maintained by only one 
point representing a Sa type galaxy. Moreover, there are no 
galaxies earlier than the Sa in the sample by Kregel et al. 
( 2002 ). 

(iii) We also demonstrate that there is no correlation be¬ 
tween the disc relative thickness and the galaxy colour. This 
fact strengthens our conclusion that there is no correlation 
between the thickness of the disc and the morphological type 
of the galaxy. 
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APPENDIX A: THE ESTIMATION OF THE 
GALAXY INCLINATION 

In this work we use two different approaches to estimate 
nearly edge-on galaxy disc inclination which are presented 
in detail by Mosenkov (2015, in prep.). 

(1) The first approach is by using the 
ExponentialDisk2>D function in the IMFIT code (Er¬ 
win 2014). The code has an advantage over GALFIT as it 
includes several 3D functions and the line-of-sight integra¬ 
tion through a 3D luminosity-density model to create a 
projected 2D image. These functions are numerically inte¬ 
grated in each image pixel. The decomposition via the 3D 
integration is much slower than the 2D edge-on disc surface 
brightness evaluation (1). To speed up computations, we 
mask the central, bulge-dominated part in galaxies and use 
the ExponentialDiskUD function in the IMFIT to fit the 
periphery of the galaxies where the disc dominates. We 


chose the radius of the masking circle to be a half of the 
radius of the galaxy, which appears to be sufficient even for 
early-type galaxies. 

The model parameters in the ExponentialDiskSD func¬ 
tion include the disc exponential scale length h, the verti¬ 
cal scale height zo, the central luminosity density 1(0,0) = 
and the inclination angle of the disc i. In order 
to find the initial estimates of the galactic disc parameters, 
we use the DECA code where the disc is described by (1). 
These estimates are assumed as the initial values for running 
IMFIT, in combination with the assumption of i = 90°. The 
output of this code is the set of the disc structural parame¬ 
ters, including the estimated inclination angle i. 

In order to verify how well this method can estimate 
the inclination angle i of the stellar disc, we performed 
the same simulations as in Section 2.2, in which instead 
of using (1) for describing the edge-on disc, we used the 
ExponentialDisk2>D function with three different values of 
i: 80°, 85° and 90°. Then we fitted these simulated images 
with the ExponentialDiskSD function of IMFIT, where the 
initial disc parameters were estimated with the DECA code, 
and the central region with the bulge was masked. The re¬ 
sults of the inclination angle estimation using IMFIT (the 
median value (i) and the standard deviation Ai) are com¬ 
pared with a priori known angles i: 

j = 80° : (i) = 81.°6, Ai = 3.°6, 

i = 85° : (i) = 85.°4, Ai = 2.°5, 

i = 90° : (i) = 88.°9, Ai = l.°8. 


The above results prove that this method can be used 
to estimate the inclination angle in the disc galaxies with¬ 
out dust attenuation and/or in the near- and mid-infrared 
bands. We found that masking of the bulge essentially re¬ 
duces the computational time while has no strong effect on 
the results of the disc decomposition. 

(2) The second method can be applied to the galaxies 
with visible dust lanes or with slight dust traces (similar to 
Bizyaev & Kajsiu (2004)). We make a natural assumption 
that the mid-planes of the stellar and dust discs coincide 
in the galaxies. We analyze ID vertical photometric profiles 
cut along the minor axis and fit them with the Sersic func¬ 
tion. The deficiency of light found by comparing the model 
and the observing profiles, and the asymmetry of the pro¬ 
file relative to the center of the galaxy indicate where the 
dust lane locates (its shift A^dust with respect to the galac¬ 
tic mid-plane). Knowing the diameter of the galaxy D and 
Azdust, we can estimate the inclination angle of the galaxy 
i (see Fig. Al) as: 


i 


arccos 


/ 2 A2:dust \ 

J ■ 


(Al) 
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Figure Al. A sketch of a nearly edge-on galaxy with the dust 
lane (see the text). 
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